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Relationships  between  the  satellite-derived  diffuse  attenuation  coefficient  of  downwelling  irradiance 
( Kd )  and  airborne-based  vertical  attenuation  of  lidar  volume  backscattering  ( a )  were  examined  in  two 
coastal  environments.  At  1.1  km  resolution  and  a  wavelength  of  532  nm,  we  found  a  greater  connection 
between  a  and  Kd  when  a  was  computed  below  2  m  depth  (Spearman  rank  correlation  coefficient  up  to 
0.96),  and  a  larger  contribution  of  Kd  to  a  with  respect  to  the  beam  attenuation  coefficient  as  estimated 
from  lidar  measurements  and  Kd  models.  Our  results  suggest  that  concurrent  passive  and  active  optical 
measurements  can  be  used  to  estimate  total  scattering  coefficient  and  backscattering  efficiency  in  waters 
without  optical  vertical  structure.  ©  2011  Optical  Society  of  America 
OCIS  codes:  010.4450,  280.0280. 


1.  Introduction 

The  diffuse  attenuation  coefficient  of  downwelling 
irradiance  (Kd)  is  a  key  optical  property  linked  to  the 
variability  of  underwater  light  fields  in  aquatic  en¬ 
vironments  [1].  For  this  reason,  Kd  has  often  been 
used  by  modelers  to  estimate  the  depth  of  the  eupho- 
tic  zone  (i.e.,  the  depth  at  which  irradiance  is  1%  of 
surface  value)  [2].  Also,  Kd  has  been  commonly  used 
to  calculate  solar  heat  budgets  [3],  determine  light 
thresholds  for  prey-predator  relationships  [4], 
estimate  coral  reef  mortality  due  to  thermal  stress 
[5],  and  indicate  water  quality  status  in  coastal 
studies  [6]. 

The  magnitude  and  vertical  distribution  of  Kd 
are  determined  by  the  sunlight  geometry  and  the 
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inherent  optical  properties  of  the  water  [7].  In  simple 
terms,  Kd  is  directly  related  to  the  total  (water+ 
particulates)  scattering  (6)  and  absorption  coefficient 
(a),  and  inversely  related  to  the  average  cosine  of 
the  zenith  angle  of  refracted  solar  photons  (direct 
beam)  just  beneath  the  sea  surface  (/*0)  (Table  1). 
Far  from  the  sea  surface,  the  Kd  distribution  is 
mainly  driven  by  variations  on  the  absorption  co¬ 
efficient  [8].  Attenuation  of  the  lidar  volume  back- 
scattering  with  depth  (a)  is  fully  or  partially  linked 
to  Kdy  depending  on  the  lidar  spot  size  at  the  sea 
surface  (R  =  H^receiver,  where  H  is  the  lidar  carrier 
altitude  above  the  sea  surface  and  0reCeiver  the 
receiver’s  field  of  view)  and  the  beam  attenuation 
coefficient  (c  =  a  +  b)  [8].  Assuming  a  vertically 
homogeneous  distribution  of  inherent  optical  proper¬ 
ties,  when  cR  «:  1,  the  exponential  decrease  of  lidar 
volume  backscattering  (S)  with  depth  is  principally 
explained  by  single  scattering  and  a  c  [9]: 
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Table  1.  List  of  Acronyms 


Symbol 

Definition 

Units 

Kd 

diffuse  attenuation  coefficient  of  downwelling  irradiance 

m1 

a 

absorption  coefficient 

m'1 

b 

scattering  coefficient 

m_1 

c 

beam  attenuation  coefficient 

m-1 

Vo 

average  cosine  of  solar  zenith  angle  beneath  the  sea  surface 

rad 

o. 

solar  zenith  angle 

rad 

a 

lidar  attenuation  coefficient 

m"1 

R 

lidar  spot  size  at  the  sea  surface 

m  rad 

H 

lidar  carrier  altitude  above  the  sea  surface 

m 

^receiver 

receiver’s  field  of  view 

rad 

s 

lidar  volume  backscattering 

m"1  sr-1 

Q 

lidar  pulse  energy 

mJ 

^rcv 

area  of  the  receiver 

m 

T 

1  atm 

transmission  of  the  atmosphere 

dimensionless 

tbw 

transmission  of  the  air/water  interface 

dimensionless 

lidar  volume  backscattering  at  180° 

m-1  sr"1 

( 

lidar  range 

m 

V 

speed  of  light  in  the  vacuum 

ms"1 

m 

refractive  index  of  seawater 

dimensionless 

bb 

backscattering  coefficient 

m"1 

bb 

backscattering  efficiency 

dimensionless 

remote  sensing  reflectance 

sr"1 

Z(001Ed(0+)) 

depth  at  which  downwelling  irradiance  is  1%  of  surface  value 

m 

Ps 

average  cosine  of  scattering 

dimensionless 

0, 

scattering  angle 

rad 

Ws) 

scattering  phase  function 

dimensionless 

S(C)  -  QArcvTltmTlviP(n)[(v/m)/(2m2(H 

+  C/m)2]exp(-2C(a  +  6)),  (1) 

where  f  is  the  lidar  range,  Q  is  the  pulse  energy,  Arcv 
is  the  area  of  the  receiver,  Tatm  and  Taw  are  the  trans¬ 
mission  of  the  atmosphere  and  the  air/water  inter¬ 
face,  respectively,  fi(jr)  is  the  volume  backscattering 
evaluated  at  a  scattering  angle  of  180°,  v  is  the  speed 
of  light  in  vacuum,  and  m  is  the  refractive  index  of 
seawater.  Conversely,  when  cR  »  1,  multiple  scatter¬ 
ing  dominates  the  received  signal,  and  Eq.  (1)  is  no 
longer  a  good  approximation  due  to  the  effects  of  vo¬ 
lume  scattering  function  shape  (i.e.,  forward  versus 
backward  directions)  and  variations  associated  with 
the  transmitter  beam  width.  In  this  case,  a  -*  Kd, 
and  the  lidar  volume  backscattering  can  be  modeled 
according  to  the  following  expression  [9]: 

S{0  *  QArcvr2tm7lw^)[(l;/m)(2m2(H  +  f/m)2] 

xexp(-2£a),  (2) 

exp(-fa)  ~  exp(-f(a  +  6))  +  exp(-af)(l  -  exp(-6f)) 
/(1  +  acr2(v/m)/ii),  (3) 

where  n  and  o2  are  the  mean  and  variance  of  the 
Gamma  probability  density  function  of  photons  as 
a  function  of  lidar  range  and  time.  The  cR  value 
corresponding  to  the  transition  between  the  “single- 
backscattering”  and  “multiple-forward-  scattering” 
regimes  is  still  in  debate  due  to  differences  between 


lidar  models  [8,9].  One  way  to  study  the  cR  threshold 
for  a  specific  lidar  system  is  to  compare  simultaneous 
a  values  with  field  measurements  of  inherent  optical 
properties  [10].  This  method  works  when  lidar  and 
optical  passive  observations  are  concurrent,  and 
have  associated  a  minimum  measurement  error  due 
to  instrument  self-shading  effects. 

Unlike  Kdl  optical  properties  influenced  by  for¬ 
ward  scattering  of  photons  (e.g.,  c,  volume  scattering 
function,  and  backscattering  efficiency  [i.e.,  bb  — 
bb/b,  where  bb  is  the  total  (water  +  particulates) 
backscattering  coefficient]  are  difficult  or  impossible 
to  study  based  on  remote  sensing  reflectance  (R^) 
signatures  [11].  For  this  reason,  most  of  studies  re¬ 
porting  6,  bb ,  and  c  distributions  in  surface  oceanic 
and  coastal  waters  rely  on  more  accurate  methods 
based  on  in-water  determinations  [12,13].  Although 
relatively  accurate  (~15%)  [14,15],  retrievals  of  Kd> 
a,  and  bb  based  on  spaceborne  or  airborne  passive 
sensors  are  not  vertically  resolved,  thus  they  re¬ 
present  integrated  values  within  the  first  optical 
depth.  This  depth  restriction  is  alleviated  when  ac¬ 
tive  optical  systems,  such  as  lidars,  are  used  instead. 
However,  Kd  values  computed  from  lidar  volume 
backscattering  profiles  are  commonly  quantified 
with  fewer  wavelengths  with  respect  to  those  Kd  va¬ 
lues  derived  from  R measurements  obtained  with 
passive  optical  systems.  Another  difference  to  em¬ 
phasize  is  that  a  is  not  necessarily  equivalent  to 
Kd  due  to  the  variable  contribution  of  forward  scat¬ 
tering  at  the  lidar  receiver.  Therefore,  for  a  specific 
light  wavelength,  a  provides  additional  information 
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relative  to  Kd  that  can  be  exploited  in  combination 
with_  passive  measurements  to  extract  c,  6, 
and  bb. 

The  aim  of  this  study  is  to  investigate  how  a  values 
derived  from  an  airborne  backscattering-based  lidar 
(i.e.,  the  National  Oceanic  and  Atmospheric  Admin¬ 
istration’s  Fish  Lidar  Oceanic  Experimental  system) 
relate  to  Kd  values  computed  from  spaceborne 
measurements  having  a  moderate  spatial  resolution 
(~1.1  km)  (Objective  1),  and  to  apply  these  measure¬ 
ments  to  estimate  b ,  bbl  and  c  within  the  first  optical 
depth  of  two  coastal  areas  (Oregon/Washington 
and  the  Afgonak/Kodiak  Shelves)  characterized  by 
waters  having  different  optical  composition  (Objec¬ 
tive  2).  Because  of  the  field  of  view  of  our  lidar  recei¬ 
ver  and  the  relatively  high  turbidity  of  the  waters 
under  investigation,  we  hypothesize  a  substantial 
correlation  between  satellite-derived  Kd  and  lidar- 
based  measurements  of  a. 

2.  Experiments 

A.  Study  Areas  and  Sampling  Design 

Comparisons  between  airborne  lidar  and  spaceborne 
passive  optical  measurements  were  performed  in  the 
eastern  shelf  of  Afgonak/Kodiak  Islands  (57.68°- 
58.03  °N,  152.04°-151.26°W)  on  17  August  2002, 
and  in  the  western  shelf  off  the  Oregon/Washington 
coast  (124.11°-124.77°W  46.16°-46.17  °N)  on  24 
August  2005  (Fig.  1).  In  the  northern  part  of  the  Gulf 
of  Alaska,  the  aerial  transect  (hereafter  AK1)  encom¬ 
passed  a  total  distance  of  40  km  and  covered  shelf 
locations  having  relatively  high  concentration  of 
chromophoric  dissolved  organic  matter  as  inferred 
from  i?rs  ratios  [16].  The  Oregon/Washington  aerial 
transect  (hereafter  OR1)  was  longer  (48  km)  and  was 
conducted  perpendicular  to  the  coast  and  along  the 
west-east  direction.  Water  properties  in  this  study 
area  are  determined  not  only  by  phytoplankton  but 
also  by  sediments  and  colored  dissolved  matter  de¬ 
rived  from  freshwater  plumes  associated  with  the 
Columbia  River.  In  general,  large  vertical  differences 
on  seawater  density  (up  to  0.22  kg/m)  reflected  a 
substantial  stratification  of  the  water  column  during 
our  aerial  surveys.  AK1  and  OR1  locations  were  se¬ 
lectively  chosen  based  on  the  following  criteria:  (1)  no 
flight  turns,  (2)  flight  track  along  the  main  orienta¬ 
tion  of  phytoplankton  or  bathymetric  gradients, 
and  (3)  locations  overlapping  satellite  pixels  with 
minimum  solar  glint.  Given  the  influence  of  flying 
altitude  on  lidar  footprint  and  signal  attenuation 
with  depth,  flight  missions  were  always  conducted 
at  a  constant  height  of  300  m. 

B.  Datasets 

1.  Active  Optical  Airborne  Measurements 

The  laser  system  used  in  this  study  was  a  nonscan¬ 
ning,  radiometric  lidar  [17]  with  three  major  compo¬ 
nents:  (1)  the  laser  and  beam-control  optics,  (2)  the 
receiver  optics  and  detector,  and  (3)  the  data  collec- 


124.5  124  0 

Longitude  (*W) 


Fig.  1.  (Color  online)  Geographic  location  of  lidar  surveys, 
(a)  Eastern  shelf  of  Afgonak/Kodiak  Islands,  (b)  Western  shelf 
off  the  Oregon/Washington  coast.  AK1  and  OKI  are  the  airborne 
transects  matching  1.1  km  satellite  ocean  color  data  (white  dots), 
land  (white),  and  bathymetry  with  1/30  deg  resolution  (color  con¬ 
tours),  missing  data  (black  areas).  Flight  direction  is  indicated 
based  on  initial  (s)  and  final  (e)  sampling  locations. 


tion  and  display  computer.  The  system  also  provides 
information  about  aircraft  position  and  attitude.  The 
laser  was  a  frequency -doubled,  Q-switched  Nd:YAG 
laser  that  produced  120  mJ  of  green  (532  nm)  light  in 
a  12  ns  pulse  at  a  rate  of  30  pulses  per  second.  The 
lidar  was  nearly  normal  to  the  sea  surface  (15°  off 
nadir),  with  a  beam  divergence  and  receiver  field  of 
view  of  17  mrad.  Signal  arriving  at  the  receiver  was 
polarized  in  the  perpendicular  direction.  Lidar  back- 
scattering  measurements  in  each  study  area  were 
collected  between  1  and  3:30  p.m.  to  maximize  tem¬ 
poral  matchup  with  spaceborne  imagery  obtained  by 
the  Moderate  Resolution  Imaging  Spectroradiometer 
(MODIS)  sensor  onboard  the  Aqua  satellite.  The 
time  difference  between  lidar  and  MODIS  measure¬ 
ments  never  exceeded  25  min. 
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2.  Passive  Optical  Spaceborne  Measurements 

MODIS-Aqua  images  (local  area  coverage,  Level  2 
OC,  1.1km  footprint)  were  obtained  from  NASA 
(http://oceancolor.gsfc.nasa.gov/).  Geo-located,  cali¬ 
brated,  and  atmospheric  corrected  R ra  values  were 
obtained  for  spectral  bands  9  (438-448  nm),  10 
(483-493  nm),  11  (526-536  nm),  12  (546-556  nm), 
and  13  (662-672  nm),  and  used  to  derive  inherent  op¬ 
tical  properties  (see  Subsection  2.C.1).  Unlike  other 
ocean  color  sensors  with  intermediate  spatial  resolu¬ 
tion  (e.g.,  SeaWiFS),  MODIS  has  a  radiometric 
channel  dedicated  to  ocean  applications  centered 
at  531  nm  (i.e.,  band  11)  and  spectrally  close  to  the 
laser  wavelength  used  in  this  study.  Also,  detection 
limits  of  MODIS  are  relatively  low  (e.g.,  signal/noise 
at  490 nm-  twofold)  with  respect  to  other  sensors. 

3.  Ancillary  Information 

Wind  speed  and  direction  during  the  AK1  and  OR1 
surveys  were  obtained  from  meteorological  stations 
located  at  the  Kodiak  Island  and  Aberdeen  airports, 
respectively  (http://www.wunderground.com/).  Recei¬ 
ver  s  radiance  contributions  due  to  foam,  glint,  and 
bubbles  are  highly  dependent  on  wind  field  charac¬ 
teristics  [181  Therefore,  wind  information  is  critical 
to  quantify  these  nonwater  radiance  contributions 
and  obtain  accurate  estimates  of  remotely  sensed 
optical  properties. 

C.  Processing  of  Remote  Sensing  Data 

1.  Satellite  Remote  Sensing  Reflectance 

The  operational  atmospheric  correction  for  ocean 
color  products  suggested  by  NASA  is  based  on  the 
Gordon  and  Wang  algorithm  [19].  The  performance 
of  this  model  may  be  compromised  in  coastal  waters 
like  those  investigated  here.  This  issue  was  exam¬ 
ined  in  our  study  areas  by  using  the  L2  OC  flag  12 
or  “TURBIDW”  MODIS-derived  total  absorption 
and  backscattering  coefficients  at  532  nm  were  com¬ 
puted  based  on  Rrs  values  at  five  wavelengths  using  a 
new  version  of  the  quasi-analytical  inversion  model 
of  Lee  et  al.  [20].  Seawater  backscattering  and  ab¬ 
sorption  coefficients  were  obtained  from  Smith  and 
Baker  tables  [21].  The  uncertainty  of  a(532)  and 
6^(532)  estimates  using  this  quasi-analytical  para¬ 
meterization  is  about  ±10%  [15].  £<*(532)  values 
were  computed  by  two  methods. 

Method  I  [22]: 

Kd( 532)  2  [a(532)  +  66(532)]/^0,  (4) 

Method  II  [14]: 

Kd( 532)  £  m0a(532)  +  mxbb{ 532) 

x  (1  -  ra2  exp(-m3a(532)).  (5) 


where  m0,  mb  ra2,  and  m3  are  parameters  that  vary 
with  solar  zenith  angle  and  depth  interval  used  to 
estimate  Kd.  Unlike  Eq.  (4),  Kd  estimates  with  Meth¬ 
od  II  are  calculated  with  a  semiempirical  parameter¬ 
ization  derived  from  radiative  transfer  theory  [14]. 

2.  Airborne  Lidar  Backscattering 

Raw  lidar  data  in  volts  were  converted  to  photo¬ 
cathode  current  values  based  on  the  specific  gain 
of  the  photomultiplier.  Afterward,  the  depth  of  each 
lidar  sample  was  found  using  the  surface  return  as  a 
reference,  and  the  return  was  multiplied  by  a  calibra¬ 
tion  factor  to  convert  photocathode  current  to  lidar 
volume  backscattering  measurements  in  m"1sr"1. 
The  calibration  factor  involves  several  parameters 
related  to  geometry  (e.g.,  sampling  altitude)  and  li¬ 
dar  system  characteristics  (e.g.,  pulse  energy,  area 
of  receiver). 

Calculation  of  a  was  performed  for  each  lidar  shot, 
followed  by  screening  and  removal  of  shots  contain¬ 
ing  subsurface  scattering  layers  [23].  The  last  step 
was  necessary  to  remove  the  influence  of  vertical 
structure  on  Kd  and  a  comparisons.  For  each  lidar 
waveform,  a  was  computed  by  linear  regression  as 
the  slope  between  water  depth  (independent  vari¬ 
able)  and  ln(S)  (dependent  variable).  This  analysis 
was  performed  using  full  resolution  profiles  (i.e., 
every  0.1  m,  accuracy  =  0.001  m)  over  different  depth 
ranges  (0-1,  0-5,  0-10,  0-15,  0-20,  2-5,  2-10,  2-15, 
and  2-20  m)  to  evaluate  the  influence  of  surface  ef¬ 
fects  (e.g.,  bubbles,  waves)  on  Kd  -  a  correlations 
and  find  the  optimum  vertical  interval  to  match 
the  penetration  depth  of  passive  sensors  (i.e.,  1  /Kd). 
Lidar  probing  depth  was  estimated  as  the  depth  at 
which  lidar  volume  backscattering  first  fell  below 
the  level  of  the  background  light  plus  10  standard  de¬ 
viations  of  the  noise  [23].  Last,  we  related  differences 
between  £<*(532)  derived  from  MODIS  at  1.1km  re¬ 
solution  with  the  median,  and  different  averages  (ar¬ 
ithmetic,  geometric,  and  harmonic)  of  a  value  within 
the  satellite  footprint.  This  numerical  exercise  was 
intended  to  examine  potential  changes  of  averaged 
a  at  1.1km  due  to  statistical  distribution  changes. 

3.  Modeling  of  b,  bb/b ,  and  c 

In  marine  waters  with  bb/a  up  to  ~0.25,  the  diffuse 
attenuation  coefficient  of  downward  irradiance  can 
be  accurately  (i.e.,  ~5%  relative  error)  approximated 
using  a,  6,  and  the  solar  altitude  [7]: 

Kd  =  /'oV2  +  G(fj0)ab)05,  (6) 

G(fd0)  =  0.425//o-0.19  for  0  S  z  £  z(o.oiEd(o+)).  (?) 

where  G  is  a  coefficient  determining  the  relative 
contribution  of  scattering  to  vertical  diffuse  attenua¬ 
tion  of  irradiance  and  is  defined  for  a  water  depth  (z) 
interval  coinciding  with  the  euphotic  zone,  i.e.,  z  at 
which  surface  downwelling  irradiance  [Ed(0±)]  is 
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reduced  in  99%.  Note  that  Eq.  (7)  (hereafter  Gl)  was 
developed  from  measurements  made  in  San  Diego 
Harbor  with  a  very  narrow  spectrum  of  6-normalized 
volume  scattering  function  or  scattering  phase  func¬ 
tion  [24].  Based  on  a  more  realistic  set  of  Monte  Carlo 
simulations  using  12  different  scattering  phase  func¬ 
tions  and  a  broader  range  of  bb/a  (0.4-2. 6)  values,  a 
new  approximation  for  G  was  found  (hereafter  G2) 
[1].  Unlike  Gl,  this  model  is  influenced  not  only  by 
geometry  of  surface  illumination  but  also  by  under¬ 
water  light  field  distribution: 


G{tto-Hs) 

=  ^0(2.636//is-  2.447) 

-(0.849//J,-  0.739), 

(8) 

i*  s  =  0.5  1 
./< 

/?(0s)cos0ssin0sd6>s, 

) 

(9) 

where  Qs  is  the  scattering  angle  in  radians  and  ps 
is  the  average  cosine  of  single  scattering  events  in 
all  directions.  In  other  words,  ps  is  a  parameter  re¬ 
lated  to  the  “shape”  of  the  volume  scattering  function 
and  can  be  empirically  linked  to  the  backscattering 
efficiency  [25]: 

jis  =  (1  -  4feft)/(1.0144  +  2.630766)  -  1.2772 b2h), 
r2  >  0.82.  0.0022  <  ft*  <0.146.  (10) 

This  relationship  is  spectrally  independent  and 
was  developed  with  869  comparisons  between  volume 
scattering  functions  and  inherent  optical  properties 
encompassing  a  broad  range  of  optical  water  types, 
and  having  b  varying  between  0.008  and  10  m”1. 
Assuming  that  a  =  Kd  and  given  that  b  =  bb/bb  with 
bb  estimated  from  inversion  modeling  [20],  we  can  cre¬ 
ate  synthetic  6* s  from  bb  values.  The  iterative  numeric 
procedure  converges  when  the  right  part  of  Eq.  (6) 
is  within  ±10%  of  a.  This  approach  was  applied  to  dif¬ 
ferent  G  functions  [Eqs.  (7)  and  (8)],  and  resulting 
6(532)  estimates  were  later  added  to  a(532)  values 
computed  from  inversion  modeling  [20]  in  order  to 
calculate  c(532). 

D.  Statistical  Analysis 

The  relationship  between  a  and  MODIS-derived 
Kd( 532)  measurements  at  1.1km  resolution  was 
quantified  using  nonparametric  Spearman  rank 
correlation  coefficients  (ps ).  The  relative  importance 
of  size  distribution  and  composition  of  particulates 
on  bb  variability  was  estimated  by  calculating  ps 
between  bb  and  two  Rrs  ratios  (R1  =  72rs(443)/ 
7?rs(488)  [26],  and  R2  =  i?rs(667)/Rrs(551))  that  are 
sensitive  to  variations  on  particle  size  distribution. 
Unlike  Rl,  R2  is  based  on  a  particle  size  distribution 
proxy  developed  with  SeaWiFS  spectral  channels 
[27].  Relationships  among  bb ,  Rl,  and  R2  were  exam¬ 
ined  using  in  situ  measurements  obtained  from  sur¬ 
face  waters  (i.e.,  0.6  m  depth)  adjacent  to  Scripps 
Institution  of  Oceanography  (University  of  California 


San  Diego).  R was  calculated  from  upwelling  radi¬ 
ance  below  the  sea  surface  and  downwelling  irradi- 
ance  above  the  sea  surface  measurements  obtained 
with  a  Hyperspectral  (wavelength  =  400-800  nm, 
spectral  resolution  =  1  nm)  Tethered  Spectral 
Radiometer  Buoy  (Satlantic  Inc.)  [28].  Particle  size 
distribution  spectra  were  characterized  with  a 
Coulter  Counter  Multisizer  III  (Beckman  Coulter, 
size  range  =  2-60  pm)  and  a  laser  diffractometer 
LISST-lOOx  (Sequoia  Scientific  Inc.,  size  range  — 

1- 200 //m).  We  quantified  the  response  of  Rl  and 
R2  as  a  function  of  the  particle  size  distribution  slope 
(Af)  estimated  from  Multisizer  III  (/)  and  LISST  (/) 
using  linear  regression  [ln(iV(D))  =  Mln(D)  ±7, 
where  N  is  the  number  of  particles  per  bin  and  unit 
of  volume  in  cubic  meters,  D  is  the  diameter  range 
in  meters,  and  I  is  the  intercept  of  the  regres¬ 
sion  curve]. 

3.  Results 

A.  Comparisons  Between  a  and  Kd 

Relative  absorption  versus  backscattering  of  pho¬ 
tons,  as  reflected  by  a(532)/6fe(532)  ratios,  varied 
between  AK1  (range,  0.040-0.052;  median,  0.046) 
and  OR1  (range,  0.035-0.091;  median,  0.057)  sur¬ 
veys;  however,  these  differences  did  not  have  a  clear 
impact  on  absolute  \a  -  Kd(532)\  values  computed  at 
1.1  km  resolution  and  based  on  a  lidar  depth  interval 
having  a  maximum  ps  between  a  and  Kd(532)  as  de¬ 
rived  from  Method  I  (Table  2).  In  general,  depth  in¬ 
tervals  with  the  highest  correlation  coefficients  were 
computed  at  depths  <5m  and  were  larger  {ps  in¬ 
creases  up  to  19.7%  in  OR1  and  14.3%  in  AK1)  when 
a  was  calculated  below  the  first  2  m  of  the  water 
column.  Also,  a  consistent  observation  at  all  depth 
intervals  under  investigation  was  the  greater  corre¬ 
lation  between  a  and  7f</(532)Method  1  in  OR1  with 
respect  to  AK1  surveys. 

The  larger  a  -  7Crf(532)Meihod  1  correlation  coeffi¬ 
cient  in  OR1  corresponded  with  a  larger  penetration 
depth  of  the  lidar  signal  and  a  larger  first  optical 
depth  (i.e.,  l/Kd\532)Meihod  l)  as  derived  from  ocean 
color  data.  In  this  area,  the  lidar  penetration  depth 
averaged  34.4  m,  based  on  a  noise  threshold  between 
5  and  15  x  10'8  sr”1  m”1.  The  penetration  depth 
varied  from  12  m  near  the  coast  to  60  m  in  those  wes¬ 
ternmost  locations  characterized  by  more  oceanic 
waters.  The  inverse  of  7C^(532)Method  1  averaged 
7.3  m,  and  varied  between  2.5  and  14.3  m  from  the 
coast  to  the  west.  The  highest  ps(a  -  7C(/(532)Method  l) 
was  obtained  over  depth  intervals  of  0-10  and 

2- 10  m  when  all  locations  (number  of  observations 
observations  -  n  —  44)  were  part  of  the  analysis. 
We  also  considered  a  reduced  dataset  (i.e.,  n  —  16) 
that  included  only  data  £30  km  from  the  coast  where 
there  was  minimal  terrigenous  material.  These  data 
had  the  greatest  correlation  at  greater  depths 
(p8  =  0.89  for  0-20  and  2-20  m).  However,  for  a  sub¬ 
sequent  interpretation  of  a  -  7Crf(532)Method  11  rela¬ 
tionships  and  the  horizontal  variability  of  a  along 
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Table  2.  Correlation  Between  Passive  and  Active  Optical  Properties  in  Oregon/Washington  and  Alaskan  Coastal  Waters* 


Experiment 

Depth  Range  (m) 

a(532) 

M532) 

ATrf(532)Method  1 

OR1 

0-1 

-0.70  (10-7) 

-0.88  (lO'7) 

-0.70  (10-7) 

0-5 

0.79  (10-7) 

0.63  (10-6) 

0.79  (lO"7) 

0-10 

0.96  (10-7) 

0.87  (lO"7) 

0.96  (lO  7) 

0-10'’ 

0.87  (10-7) 

0.90  (lO'7) 

0.87  <10-7) 

0-15'1 

0.87  (10-7) 

0.87  (10-7) 

0.87  (lO"7) 

0-20'’ 

0.89  (10-7) 

0.90  (10-7) 

0.89  (10-7) 

2-5 

0.94  (10-7) 

0.83  (lO’7) 

0.94  (10-7) 

2-10 

0.96  (10-7) 

0.87  (10‘7) 

0.96  (lO  7) 

2-10* 

0.86  (10-7) 

0.88  (10-7) 

0.86  (10-7) 

2-15* 

0.87  (10-7) 

0.87  (lO’7) 

0.87  (10-7) 

2-20* 

0.89  (10-7) 

0.90  (10-7) 

0.89  (10-7) 

AK1 

0-1 

-0.64  (9 10-5) 

-0.673  (4  lO'5) 

-0.662  (2  lO'5) 

0-5 

2.8 10-3  (0.99 )"” 

0.025  (0.859)ns 

-3.6  lO’3  (0.998)ns 

0-10 

0.64  (210-5) 

0.700  (2  lO’5) 

0.642  (3 10-5) 

0-15 

0.74  (2  lO'5) 

0.745  (2 10-5) 

0.736  (210-7) 

0-20 

0.70(2  lO'5) 

0.706  (2  lO  3) 

0.710  (210-7) 

2-5 

0.64  (2  lO’5) 

0.664(6  lO'6) 

0.649  (2  lO"5) 

2-10 

0.73  (2  lO'7) 

0.771  (2 10-7) 

0.734  (2  lO'7) 

2-15 

0.76  (210'7) 

0.756  (2 10-7) 

0.757  (2  lO'7) 

2-20 

0.70  (2  lO'7) 

0.703  (2  lO'7) 

0.704  (2  lO"7) 

“For  each  correlation,  probability  of  accepting  the  null  hypothesis  ( Hl} ,  p„  =  0;  i.e.,  variables  are  uncorrelated)  is  indicated  between 
parentheses.  Nonsignificant  correlations  at  95%  (ns)  confidence  level,  highest  a  -  /frf(532)Method  1  correlations  are  highlighted  in  bold. 
OR1  and  AK1  are  defined  in  Subsection  2. A. 

^Calculated  with  measurements  obtained  at  £30  km  distance  with  respect  to  the  starting  flying  point. 


each  transect,  we  used  a  calculated  between  2  and 
10  m  to  provide  more  comparisons  and  representa¬ 
tive  samples  of  the  whole  survey.  In  AK1,  the  pene¬ 
tration  depth  of  the  lidar  averaged  20  m  for  a  noise 
threshold  between  3  and  11  x  10~10  sr_1  m"1.  The  op¬ 
tical  depth  was  generally  shallower  ( x  =  5.9  m)  and 
less  variable  (4. 1-7.5  m)  than  at  OR1.  The  maximum 
/?s(cr-^(532)Meth(KiI1)  was  obtained  for  a  values 
computed  using  a  2-15  m  depth  interval  (p8  =  0.757, 
n  =  36). 

The  similarity  between  Kd  (532 )  and  a  computed  be¬ 
low  2  m  depth  was  corroborated  based  on  comparisons 
between  a  and  K Method  11  (Table  3).  Since  the  Method  II 
algorithm  only  resolves  Kd(532)  at  three  sunlight  geo¬ 
metries,  /^Method  II  values  were  calculated  between 
two  solar  zenith  angles  (i.e.,  60°  and  30°)  in  order 
to  include  all  possible  solar  elevations  during  the  air¬ 
borne  campaigns  (i.e.,  42°  to  48°).  As  expected,  abso¬ 
lute  differences  between  the  magnitude  of  a  computed 
at  the  maximum  ps(a-Kd( 532)MethodI)  and  Kd 
(532)Method  11  tended  to  be  minimum  for  larger  depth 
ranges.  In  OR1,  the  satellite  overpass  on  24  August 
2005  was  relatively  early  (i.e.,  local  time  =  13:05  h) 
and  the  root  mean  square  error  (RMS)  between  lidar 
and  passive  attenuation  coefficients  was  smaller 
when  the  Sun  was  close  to  the  vertical  (i.e.,  Method 
I  simulations  at  30°).  Conversely,  MODIS-Aqua  i?re 
measurements  during  lidar  surveys  of  17  August 
2002  were  obtained  late  in  the  afternoon  (local  time  = 
15:05  h),  resulting  in  a  closer  agreement  between  a 
and  /fd(532)Metho°  11  at  the  two  solar  positions  (e.g., 
RMS  =  0.022,  depth  range  =  0-20  m,  zenith  angle 
30°  and  60°). 

Overall,  a  values  at  the  optimum  correlation 
depth  range  were  closer  to  /fd(532)Method  11  than 


if£/(532)Method  l.  The  relative  biases  {i.e.,  [(a-Kd 
(532))/iQ(532)])  as  percentages  were  +9.3% 
(Method  II)  and  +11.3%  (Method  I)  in  OR1,  and 
+8.7%  (Method  II)  and  +8.9%  (Method  I)  in  AK1. 
These  biases  were  much  greater  than  Kd  differences 
(<2%)  associated  to  changes  on  solar  altitude  during 
each  aerial  survey  (~0.03%).  Also,  statistical  varia¬ 
tions  of  a  within  the  satellite  footprint  had  a  minor 
impact  on  observed  Kd( 532)  -  a  differences  since  the 
magnitude  of  different  a  averages  (i.e.,  arithmetic, 
geometric,  harmonic)  and  median  only  differ  in  the 
fourth  decimal  unit. 

The  light  interaction  mechanism  explaining  the 
aforementioned  a- A'(y(532)Mothod  1  correlations  dif¬ 
fered  between  study  areas.  In  OR1,  the  lidar  attenua- 


Table  3.  Difference  Between  a  and  AC</(532)M*,hod  ■  for  Two  Different 
Solar  Zenith  Angles* 


Experiment 

Depth 
Range (m) 

0Z  -  0.5235 

ez  =  1.0471 

0-1 

0.040(17.9) 

0.048(12.2) 

OR1 

0-5 

0.038  (15.5) 

0.049  (12.4) 

0-10 

0.037  (13.5) 

0.049  (12.4) 

0-20* 

0.008  (12.0) 

0.012  (6.8) 

AK1 

0-1 

0.029  (16.3) 

0.021  (9.1) 

0-5 

0.025  (13.8) 

0.022  (9.4) 

0-10 

0.022  (11.3) 

0.023  (9.5) 

0-20 

0.022  (11.4) 

0.022  (9.3) 

“Each  value  corresponds  to  root  mean  square  between  the  arith¬ 
metic  average  of  a  at  the  respective  depth  interval  (i.e.,  OR1, 
2-10  m;  AK1,  2-15  m)  and  MODIS-derived  Kd( 532)  computed  at 
each  1.1km  pixel.  Between  parentheses  is  the  relative  difference 
as  percentage;  i.e.,  100  ((a- ^(532)Method  u)/tfd(532)Method  “]. 

*Idem  to  Table  2. 


20  June  2011  /  Vol.  50,  No.  18  /  APPLIED  OPTICS 


2995 


tion  below  2  m  depth  was  primarily  driven  by  a  (532) 
changes.  Conversely,  in  AK1,  the  effects  of  a  (532)  and 
6^(532)  on  a  were  comparable.  This  finding  is  better 
illustrated  in  Fig.  2,  where  we  evaluated  the  horizon¬ 
tal  coherence  between  satellite-based  /frf(532)Method  1 
and  airborne-based  a  for  spatial  scales  between  1  and 
50  km.  For  this  analysis,  we  used  Method  I  rather 
than  Method  II  to  derive  Kd( 532)  because  Method 
I  depends  on  //0,  while  Method  II  estimates  are  con¬ 
strained  by  solar  position. 

In  OR1,  i^(532)^ethod  x,  a(532),  and  a  decreased  by 
80%  from  the  coast  to  offshore  locations  [Fig.  2(a)]. 
On  the  other  hand,  6^(532)  was  relatively  low 
(<0.015  m"1)  within  the  first  6  km  of  the  lidar  survey, 
reached  a  maximum  (up  to  0.018  m"1)  at  intermedi¬ 
ate  distances  from  shore  (9  to  21  km),  and  decreased 
to  minimum  values  (<0.005  m"1)  by  the  offshore  end 
of  the  transect.  Not  surprisingly,  satellite-derived 
iQ(532)Method  1  estimates  had  a  larger  uncertainty 
(±20%)  than  averaged  a  values  (±5%  in  average,  n 
per  pixel  =  200)  within  the  same  1.1km  footprint. 
Despite  these  similarities,  there  were  some  sections 
along  the  transect  (e.g.,  24  to  29  km)  where 
jffd(532)Method  1  error  bars  did  not  overlap  the  arith¬ 
metic  mean  of  lidar  attenuation  coefficient  values. 
Like  the  OR1  results,  the  spatial  patterns  of 
/fc/(532)MGthod  1  and  a  in  AK1  had  a  positive  covaria¬ 
tion,  and  a  was  always  greater  than  /Q(532)Method  1 
[Fig.  2(b)].  However,  there  were  striking  differences 
in  terms  of  how  the  inherent  optical  properties 
affected  Kd.  Indeed,  6^(532)  values  in  the  Alaskan 
Shelf  covaried  with  ( 532 )Method  *,  a(532)  and  a, 
and  were  related  to  variations  in  the  bathymetry 
(i.e.,  50-70  m  depth  over  shallow  “banks”  versus 
150-250  m  over  the  “troughs”).  The  presence  of 
“banks”  (e.g.,  18-22  km  from  the  start)  and  “troughs” 
(e.g.,  10-12  km)  were  consistently  discriminated 
by  MODIS-Aqua  and  lidar  measurements  as  areas 
characterized  by  relatively  high  (1/a  =  5. 7-6.0  m, 
l/tfrf(532)Method  1  -  6.0-7. 7  m)  and  low  (1/a  = 
4. 2-4.5  m,  l//Q(532)Method  1  =  4.2-4.8  m)  water 
visibility  values,  respectively. 

B.  Analysis  of  Inherent  Optical  Properties  Computed  from 
a  and  Kd 

Assuming  that  oris  not  different  from  Kd,  we  can  use 
Eq.  (6)  to  model  bh  and  subsequently  b  and  c  based  on 


Fig.  2.  (Color  online)  Spatial  coherence  between  satellite-derived 
/£rf(532)Method  1  and  airborne-based  a.  (a)  OR1;  (b)  AK1.  For  each 
comparison,  a  is  the  arithmetic  mean  within  the  satellite  footprint 
(pink  solid  curve,  left  axis),  a (532)  (blue  solid  curve,  left  axis), 
JQ(532)Method  1  (gray  solid  curve,  left  axis),  and  66(532)  (red  solid 
curve,  right  axis),  missing  lidar  data  (hatched  bars).  To  better  il¬ 
lustrate  /frf(532)Method  1  -a  differences,  error  bars  of  a(532)  and 
6/,(532)  are  not  shown. 

active  and  passive  optical  properties.  The  approxi¬ 
mation  a  =  Kd  was  met  by  selecting  a  subset  of  va¬ 
lues  in  each  transect  having  the  minimum  \a-Kd\ 
magnitude.  Under  this  premise,  we  extracted  four 
groups  from  Fig.  2:  “coastal”  (locations  5,  9,  and 
11km)  and  “oceanic”  (32,  44,  and  46  km)  in  OR1, 
and  “banks”  (2,  4,  and  22  km  )  and  “troughs”  (9,  10, 
and  30km)  in  AK1.  Errors  in  calculating  6(532), 
6^(532),  andc(532)  from  Eq.  (6)  with  G1  and  G2  were 
24%,  26%,  38%,  and  77%,  79%,  83%,  respectively.  In 
general,  higher  6(532)  (up  to  0.81  m"1)  and  c(532)  (up 
to  1.09m"1),  and  lower  6^(532)  (up  to  0.015)  values 


Table  4.  Summary  of  Inherent  Optical  Properties  Estimated  from  a  and  Eqs.  (6>— (8)* 


Experiment 

Gl 

G2 

5(532) 

*6(5  32) 

c(532) 

6(532) 

*6<532) 

c(532) 

OR1 

Coastal 

Min 

0.494 

0.020 

0.747 

0.272 

0.015 

0.524 

Max 

0.607 

0.033 

0.877 

0.813 

0.060 

1.089 

Oceanic 

Min 

0.016 

0.031 

0.074 

0.020 

0.040 

0.077 

Max 

0.173 

0.120 

0.265 

0.074 

0.100 

0.166 

AK1 

Banks 

Min 

0.236 

0.023 

0.367 

0.040 

0.002 

0.17 

Max 

0.320 

0.025 

0.469 

4.091 

0.146 

3.42 

Troughs 

Min 

0.120 

0.018 

0.256 

0.026 

0.002 

0.19 

Max 

0.466 

0.024 

0.645 

3.273 

0.120 

4.37 

"Range  of  values  for  each  estimate  based  on  one  (i.e.,  Gl)  or  many  (i.e.,  G2)  water  types. 
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Fig.  3.  Response  of  remote  sensing  reflectance  ratios  to  variabil¬ 
ity  of  particle  size  distribution.  Rl,  /?„(443)/?„(488)  (left  axis) 
and  R2,  /?rg(670)/i?rs(555)  (right  axis)  are  plotted  as  functions  of 
particle  size  spectrum  slope  derived  from  Coulter  counter 
(a)  Multisizer  III  (/)  and  (b)  LISST  (y). 

were  estimated  in  OR1  than  in  AK1  (Table  4).  In 
OR1,  the  “coastal”  group  was  characterized  by  higher 
6(532)  (up  to  tenfold)  and  c(532)  (up  to  sevenfold) 
and  lower  6^(532)  (up  to  0.4-fold)  values  than  those 
for  the  “oceanic”  group.  Shallower  waters  of  AK1  as¬ 
sociated  with  submarine  banks  had  typically  lower 
6(532)  (up  to  0.24  m-1)  and  c(532)  (up  to  0.37  m-1) 
and  higher  6/,(532)  (up  to  0.025)  values  than  those 
over  deep  canyons  (i.e.,  “troughs”). 

The  concurrent  use  of  active  and  optical  passive 
measurements  in  this  study  was  also  applied  to 
estimate  second-order  optical  attributes  affecting 


Table  5.  Statistical  Relationships  Between  Particle  Size  Distribution 
and  Rn  Ratio  Variability® 


M 

/ 

n 

r2 

Rl 

r 

0.050  (0.018)* 

0.518  (0.048)* 

5 

0.712 

X 

0.080  (0.030)* 

0.403  (0.091)* 

5 

0.704 

R2 

r 

0.063  (0.024)* 

-0.048  (0.062) 

5 

0.692 

X 

0.073  (0.056) 

-0.108  (0.172) 

5 

0.358 

°The  linear  model  used  to  estimate  slope  of  particle  concentra¬ 
tion  (y )  as  a  function  of  particle  size  range  (x)  isy  —  Mx  +  /;  Af  and  I 
are  the  slope  and  intercept  of  the  regression  curve,  respectively. 
Between  parentheses  is  the  standard  error  of  each  regression  coef¬ 
ficient,  Af  and  /  are  different  from  0  at  95%  confidence  level  (*),  and 
n  is  the  number  of  comparisons.  Rl,  R2,  y,  and  /  are  explained  in 
Subsection  2.D. 


Table  6.  Influence  of  Particle  Size  Distribution  on  Spatial 
Variability  of  6b(532)* 


Experiment 

P* 

P 

n 

OR1 

0.60 

0.24 

6 

0.60 

0.24 

6 

AK1 

0.78 

0.06 

6 

0.87 

0.03* 

6 

aps  is  the  Spearman  correlation  coefficient,  P  is  the  probability 
of  rejecting  the  null  hypothesis  ( H0 ,  ps  —  0)  at  95%  confidence 
level  (*),  and  n  is  the  number  of  comparisons.  For  each  subset,  first 
and  second  row  correspond  to  bb  estimates  using  G1  and  G2, 
respectively. 

magnitude  of  6^(532)  in  each  marine  environment. 
Preliminary  results  based  on  data  collected  in  near¬ 
shore  waters  of  the  Southern  California  Bight 
revealed  a  positive  covariation  between  the  slope 
of  particle  size  distribution  estimated  from  two  in¬ 
struments  (Multisizer  III  and  LISST),  and  two  differ¬ 
ent  Rrs  ratios  in  the  visible  range  (Fig.  3).  Based  on 
linear  regression,  we  found  that  remote  sensing  re¬ 
flectance  indices  of  particle  size  distribution  had  a 
greater  covariation  with  particle  size  spectra  mea¬ 
sured  by  Multisizer  III.  This  relationship  was  stron¬ 
ger  (i.e.,  higher  coefficient  of  determination,  r2)  for 
comparisons  based  on  Rrs(443)/7?rs(488)  (Table  5). 
The  lower  performance  of  R2  for  detecting  changes 
on  particle  size  distribution  was  even  worse  (i.e., 
r2  up  to  14%  lower)  when  R2  was  calculated  with  Sea- 
WiFS  spectral  channels  [i.e.,  Rrs(670)/7?rs(555)]  [27]. 
Thus,  we  only  examined  a  posteriori  correlations  be¬ 
tween  7?rs(443)/i?rs(488)  and  6^(532)  values.  These 
simple  correlations  illustrate  the  greater  effect  of 
particle  size  distribution  changes  on  6*,  (532)  varia¬ 
bility  in  AK1  than  in  OR1,  and  the  greater  impact 
of  different  particle  size  ranges  on  G2-derived  inher¬ 
ent  optical  properties  than  those  derived  from  G1  in 
AK1  (Table  6). 

4.  Discussion 

The  value  of  a  in  oceanographic  studies  will  vary  due 
to  the  type  of  lidar  system  and  the  turbidity  of  the 
water  under  study.  At  1.1km  resolution,  concurrent 
airborne-derived  a  and  satellite-derived  Kd  in  shelf 
waters  of  Alaska  and  Oregon/Washington  during 
summer  showed  that  a  during  these  experiments 
was  a  good  predictor  of  Kd  in  the  green  spectral  range 
(i.e.,  wavelength  =  532  nm).  Indeed,  and  based  on 
single  linear  regression  models,  we  found  for  each 
survey  that  a  (dependent  variable)  was  related  to 
Kd  (independent  variable)  with  a  regression  slope 
(AK1,  1.074  ±0.132  standard  error,  n  —  36;  OR1, 
0.975  ±  0.053,  n  =  44)  and  intercept  (AK1,  -0.030  ± 
0.025  standard  error;  OR1,  -0.013  ±  0.013)  not  sta¬ 
tistically  different  from  1  and  0,  respectively.  The 
connection  between  lidar-derived  a  and  Kd  was 
previously  reported  in  the  Southern  California 
Bight  [10]. 

Similar  to  our  study,  they  found  that  a  values  at 
20  m  depth  were  mainly  determined  by  Kd  and  not 
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c  (see  linear  regression  slope  of  Eqs.  5  and  6  in  [10]). 
However,  in  contrast  with  our  findings,  their  a  -  Kd 
regression  slope  was  below  unity.  We  attribute  this 
apparent  discrepancy  to  differences  in  cR.  In  the 
Southern  California  Bight  study  [10],  the  laser  beam 
divergence  angle  (43  mrad)  and  receiver  field  of  view 
(26  mrad)  were  larger  with  respect  to  our  study,  and 
lidar  measurements  were  obtained  from  the  ship 
deck  (i.e.,  distance  between  lidar  source/detector  and 
sea  surface  was  10.3  m).  Given  this  geometry,  their  R 
(-0.08  m)  and  cR  (-0. 01)  values  were  relatively  small 
compared  with  our  values  (R  -  2.5  m,  cR  -  1).  There¬ 
fore,  as  cR  becomes  smaller  than  1  (i.e.,  Churnside 
et  al.'s  study  [10]),  c  is  expected  to  explain  a  larger 
fraction  of  Kd  [29].  Another  variable  decreasing  cR 
in  contribution  [10]  was  their  lower  c(532)  values 
(mean  =  0.098  m'1)  with  respect  to  our  study  (AK1, 
0.173  m'1;  OR1,  0.193  m'1).  Note  that  Kd  and  c  deter¬ 
minations  by  [10]  were  more  accurate  than  ours  since 
they  were  derived  from  in-water  measurements. 
Finally,  it  is  worth  emphasizing  that  our  lidar  results 
are  based  on  cross-polarized  lidar  returns,  while  the 
study  in  [10]  describes  values  for  the  copolarized  re¬ 
turns.  Multiple  forward  scattering  can  produce  a  re¬ 
duced  attenuation  of  the  cross-polarized  returns 
under  some  conditions  [30]. 

Despite  the  overall  agreement  between  magnitude 
of  a  and  MODIS-derived  Kd( 532)  measurements 
during  the  AK1  and  OR1  surveys,  we  detected  sub¬ 
stantial  changes  (up  to  -f 0.088  in  OR1  and  + 0.049  in 
AK1)  between  these  two  properties  at  specific  loca¬ 
tions  along  the  transects  (e.g.,  26-28  km  in  OR1, 
18  km  in  AK1).  Since  a  values  lie  between  c  and 
Kd  [8-10]  and  c  is  always  larger  than  Kdy  it  is  sug¬ 
gested  that  the  observed  increase  of  a  with  respect 
to  the  Kd( 532)  magnitude  was  associated  in  these  lo¬ 
cations  with  a  larger  relative  contribution  of  c(532) 
to  a,  and,  consequently,  a  greater  proportion  of  for¬ 
ward  scattering  defining  the  underwater  light  field. 
We  attribute  these  spatial  changes  (i.e.,  within  and 
between  transects)  to  the  presence  of  different  opti¬ 
cal  water  types. 

Although  averaged  hourly  wind  speed  was  higher 
during  the  OR1  survey  (4.44  ms"1)  than  during  the 
AK1  survey  (3.33  m  s"1 ),  its  influence  on  a  -  Kd( 532) 
differences  was  secondary  due  to  three  main  reasons. 
First,  and  based  on  other  studies  [31],  the  greater 
wind  intensity  and  associated  production  of  subsur¬ 
face  bubbles  in  Oregon/Washington  is  expected  to 
have  a  minor  impact  on  lidar  volume  backscattering 
(-25%)  compared  to  observed  relative  changes 
between  a  and  Kd( 532)  (up  to  60%).  Second,  our  a  - 
Kd{ 532)  comparisons  were  based  on  a  calculated  be¬ 
low  2  m  depth,  thus  eliminating  major  interference 
due  to  “surface”  effects.  This  interference  was  more 
pronounced  when  Kd{ 532)  was  estimated  with 
Method  I  since  the  Method  II  approach  was  devel¬ 
oped  using  constant  and  relatively  weak  winds 
(i.e.,  5  ms'1).  Last,  the  influence  of  wind-mediated 
changes  on  sea  surface  slopes  (cross  and  up/down) 
and  subsequent  contribution  of  Sun  glint  to  i?re 


was  ruled  out  as  a  possible  major  bias  of  Kd{ 532) 
estimates  since  the  specific  radiance  threshold  at 
865  nm  proposed  for  MODIS-Aqua  (i.e.,  flag 
MODGLINT  or  moderate  Sun  glint  contamination) 
was  never  exceeded  during  the  analysis  of  ocean 
color  data. 

For  a  specific  range  of  oceanographic  conditions 
and  lidar  system  parameters,  active  and  passive  op¬ 
tical  measurements  were  combined  in  this  study  to 
calculate  inherent  optical  properties  related  to  mag¬ 
nitude  and  angular  distribution  of  light  scattering. 
Values  of  6(532),  c(532),  and  6^(532)  can  also  be  es¬ 
timated  by  making  equal  Eqs.  (5)  and  (6)  and  solving 
for  6.  However,  this  mathematical  procedure  has 
some  drawbacks.  First,  the  use  of  two  models  is  add¬ 
ing  more  error  (-12%)  to  the  final  estimates.  Second, 
Kd  models  based  on  passive  optical  measurements 
are  based  on  b^/a  changes,  thus  forward  scattering 
contributions  are  not  quantified.  The  median  of 
6(532)  (n  =  3)  in  OR1  (Gl,  0.49  m'1)  and  AK1  (Gl, 
0.36  m'1)  as  estimated  from  Gl  was  within  the  range 
of  oceanic  (0.275  m-1)  and  coastal  (1.21  to  1.82  m"1) 
values  reported  at  the  same  wavelength  in  Southern 
California  waters  [24].  Our  Gl-based  c(532)  esti¬ 
mates  (median  in  OR1, 0.55  nr1;  AKl.  0.48  m"1)  were 
also  intermediate  between  oligotrophic  (e.g.,  0.2  m"1) 
and  eutrophic  (e.g.,  up  to  10  m"1)  marine  environ¬ 
ments  [32].  Unlike  6(532)  and  c(532),  our  estimates 
of  backscattering  probability  at  532  nm  were  some¬ 
times  beyond  (e.g.,  oceanic  in  OR1,  G2-based  in 
AKl)  the  maximum  bb  values  measured  in  the  Paci¬ 
fic  Central  Gyre  (0.04-0.06)  [13]  and  the  Bahamas 
Shelf  during  a  “whiting  event”  (-0.05)  [12].  Given 
the  relative  low  signal  in  our  oceanic  locations  and 
the  high  uncertainty  (relative  error  up  to  79%)  of 
6^(532)  estimates  using  more  realistic  volume  scat¬ 
tering  functions  (i.e.,  G2),  we  suggest  that  observed 
variations  of  6^(532)  with  respect  to  literature  values 
were  apparent.  The  use  of  multiple  and  different  op¬ 
tical  signals  in  this  research  allowed  us  not  only  to 
quantify  a  budget  of  inherent  optical  properties 
but  also  to  investigate  additional  physical  aspects  re¬ 
lated  to  bb  behavior  due  to  variations  on  particle  size 
distribution.  In  that  regard,  the  substantial  and  ex¬ 
clusive  correlation  found  between  G2-based  bb  and 
/?rs(443)//?rs(488)  in  AKl  is  indirectly  suggesting 
two  important  facts:  the  more  heterogeneous  nature 
of  the  volume  scattering  function  in  AKl  with  respect 
to  OR1,  and  the  greater  importance  of  other  factors 
(e.g.,  changes  on  particle  composition  due  to  river 
sediments)  affecting  backscattering  efficiency  during 
OR1  surveys. 

This  work  was  supported  by  the  Naval  Research 
Laboratory  (NRL)  internal  project  “3D  Remote 
Sensing  with  a  Multiple-Band  Active  and  Passive 
System:  Theoretical  Basis,”  PBE0601153N. 
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